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ABSTRACT 

Context. At this time there does not exist a robust set of rules connecting low and high /? waves across the y6 a 1 layer. The work 
here contributes specifically to what happens when a low fj fast wave crosses the yS « 1 layer and transforms into high j8 fast and slow 
waves. 

Aims. The nature of fast and slow magnetoacoustic waves is investigated in a finite f3 plasma in the neighbourhood of a two- 
dimensional null point. 

Methods. The linearised equations are solved in both polar and cartesian forms with a two-step Lax-Wendroff numerical scheme. 
Analytical work (e.g. small /3 expansion and WKB approximation) also complement the work. 

Results. It is found that when a finite gas pressure is included in magnetic equilibrium containing an X-type null point, a fast wave is 
attracted towards the null by a refraction effect and that a slow wave is generated as the wave crosses the /? a; 1 layer. Current accu- 
mulation occurs close to the null and along nearby separatrices. The fast wave can now pass through the origin due to the non-zero 
sound speed, an effect not previously seen in related papers but clear seen for larger values of /?. Some of the energy can now leave 
the region of the null point and there is again generation of a slow wave component (we find that the fraction of the incident wave 
converted to a slow wave is proportional to /?). We conclude that there are two competing phenomena; the refraction effect (due to the 
variable Alfven speed) and the contribution from the non-zero sound speed. 

Conclusions. These experiments illustrate the importance of the magnetic topology and of the location of the /3 « 1 layer in the 
system. 
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1. Introduction 

MHD wave motions have recently been observed in the Sun's at- 
[ mosphere with the SOHO and TRACE satellites, see for example 
, the detection of slow magnetoacoustic waves by Berghmans & 
Clette, (1999) and De Moortel et al. (2000) and of fast magne- 
toacoustic waves by Nakariakov et al. (1999). It is clear that the 
local coronal magnetic field plays a key role in determining their 
propagation properties. 

To begin to understand this inhomogeneous magnetised en- 
vironment, it is useful to look at the structure (topology) of 
\ the magnetic field itself. Potential field extrapolations of the 
coronal magnetic field can be made from photospheric mag- 
netograms. Such extrapolations show the existence of an im- 
portant feature of the topology; null points. Null points are 
points in the field where the Alfven speed is zero. Detailed in- 
vestigations of the coronal magnetic field, using such potential 
field calculations, can be found in |Beveridge et al. (2002) and 
[Brown & Priest (200 1)| 

Building on the earlier work of Bulanov and Syrovatskii (1980) 
and Craig and Watson (1992), McLaugMin & Hood (2004) 
found that for a single 2D null point, the fast magnetoacoustic 
wave was attracted to the null and the wave energy accumulated 
there. In addition, they found that the Alfven wave energy accu- 
mulated along the separatrices, the topological feature that sepa- 
rates regions with diff'erent magnetic flux connectivity. Their pa- 
per looked at MHD wave propagation in ayS = plasma. The aim 
of this paper is to extend their model to include plasma pressure 



(finite /? plasma). The most obvious efifect is the introduction of 
slow magnetoacoustic waves. The fast wave can now also pass 
through the null point (as we have a non-zero fast wave speed 
there due to the finite sound speed) and, thus, perhaps carry wave 
energy away from that area. There could also be coupling and 
wave conversion near the location where the sound speed and 
Alfven speed become comparable in magnitude. However, the 
exact nature of such coupling in MHD is unknown and will be 
looked at here. The behaviour of the Alfven waves is unaffected 
by a finite /3 (the plasma pressure plays no role in its propagation) 
and so the description by McLaughUn & Hood (2004), remains 
valid in this linear, 2D regime. 

Waves in the neighbourhood of a single 2D null point have been 



investigated by various authors. [Bulanov & Sy rovatskii (1980) 
provided a detailed discussion of the propagation of 
fast and Alfven waves using cylindrical symmetry. 
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Craig & Watson (1992)| mainly consider the radial propagation 
of the m - mode (where m is the azimuthal wavenumber) 
using a mixture of analytical and numerical solutions. They 
show that the propagation of the m - Q wave towards the null 
point generates an exponentially large increase in the current 
density and that magnetic resistivity dissipates this current in 
a time related to log//. Craig and McClymont (1991, 1993) 
investigate the normal mode solutions for both w = and m + 
modes with resistivity included. Again, they emphasise that the 
cuiTent builds up as the inverse square of the radial distance 
from the null point. All these investigations were caiTied out 
using cylindrical models in which the generated waves encircled 
the null point. In a sense, there is nowhere else for the wave to 
propagate except into the null point. 
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Fig. 1. The equilibrium equilibrium magnetic field and coordi- 
nate system. 



A very detailed and comprehensive set of 2D numerical 
simulations of wave propagation in a stratified magneto- 
atmosphere was conducted by [Rosenthal et al. (2002)| and 
Bogdan et al. (2003) In these simulations, an oscillating piston 
generated both fast and slow MHD waves on a lower boundary 
and sent these waves up into the stratified, magnetised plasma. 
Their calculations showed that there was coupling between the 
fast and slow waves, and that this coupling was confined to a thin 
layer where the sound speed and the Alfven velocity are compa- 
rable in magnitude, i.e. where the plasma-beta approaches unity. 
Away from this conversion zone, the waves were decoupled as 
either the magnetic pressure or plasma pressure dominated. In 
this, their papers and ours have a similar goal; to see how the 
topology affects the propagation of MHD waves, in a 2D sys- 
tem where the ratio of the sound speed to the Alfven speed varys 
along every magnetic line of force. 

Other authors have also looked at MHD mode coupling, 
describes 



Cally & Bogdan (1997) describes 2D simulations in which 
both /-modes and /5-modes are (partially) converted to slow 
magnetoacoustic gravity waves, due to strong gravitational 
stratification. |De Moortel et al. (2 004) investigated driving slow 
waves on the boundary of a 2D geometry with a horizon- 
tal density variation. They found coupling between slow and 
fast waves and phase mixing of the slow waves. The cou- 
pling of different wave modes has also been investigated by 
Ferraro & Plumpton (19 58), Zhugzhd & Dzha hlov (1982)| (with 
Meijer G-functions) and |Cally (2001) (with hypergeometric 2F3 
functions). All these works considered mode coupling through a 
gravitational stratification (vertical density inhomogenity) and, 
in particular, investigated the propagation of waves from a high 
jS to a low jS plasma. As waves propagate towards a null point, 
it is the nature of low p to high j6 propagation that is important. 
Finally, the coupling of fast waves and Alfven waves has been 
investigated by Parker (1991) (linear with a density gradient) and 
by Nakariakov et al. (1997), (nonlinear excitation). 
The paper has the following outline. In Section|2] the basic equa- 
tions are described and the importance of the plasma p discussed. 
The results for the fast and slow magnetoacoustic waves are pre- 
sented for a reference value of the plasma namely ySo = 0.25 
in Section |3] Section |4]provides an interpretation of our results, 
with a discussion of mode conversion in Section |5] The robust- 
ness of our results and a clearer demonstration of mode con- 
version is shown in Section|6]fory6o = 2.25. The conclusions are 
given in Section 7 and the appendices provide analytical approx- 
imations. 



2. Basic Equations 

The usual MHD equations appropriate to the solar corona are 
used, with pressure and resistivity included. Hence, 



P 



-H (V ■ V) V 



= -(VxB)xB- , 



5B 

'dt 



— = V X (v X B) -H yyV^B 



= 



dp 
dt 



+ (\ ■V)p - -y/?V ■ V , 



where p is the mass density, v is the plasma velocity, B the 
magnetic induction (usually called the magnetic field), p is the 
plasma pressure, fi - Anx 10" ■'Hm^' the magnetic permeability, 
rj - 1 //icr is the magnetic diffusivity (m's"' j and cr the electrical 
conductivity. 

2. 1 . Basic equilibrium 

The basic potential magnetic field structure is taken as a simple 
two dimensional, X-type neutral point. There are a lot of compli- 
cated effects including mode conversion and coupling; and a 2D 
geometry allows one to understand these effects better, before 
the extension to 3D. Therefore, the equilibrium magnetic field 
(see Figure [TJ is taken as 



Bo = fip,0,-£), 
Vfl a/ 



(1) 



where B is a characteristic field strength and a is the length scale 
for magnetic field variations. Obviously, this configuration is no 
longer valid far from the null point, as the field strength tends to 



infinity. However, McLaughlin & Hood (2006) looked at a mag 



netic field with a field strength that decays far from the null 
(again for a j8 = plasma) and they found that the key results 
from McLaughlin & Hood (2004)lre main valid close to the null. 
As in [McLaughlin & Hood (2004)| the Unearised MHD equa- 
tions are used to study the nature of the wave propagation near 
the null point. Using subscript for equilibrium quantities and 
1 for perturbed quantities, the Unearised equation of motion be- 



Po 



d\i /VxBi 



X Bo - V/?i 



dt \ n 
the linearised induction equation 

^ = V X (vi X Bo) + T/V^Bi , 

dt 

the linearised equation of mass continuity 

dpi 
dt 



+ V ■ (povi) = , 



(2) 



(3) 



(4) 



(5) 



and the adiabatic energy equation 

= -ypQ (V ■ vi) . 

dt 

We will not discuss equation ((H further as it can be solved 
once we know vi. We assume the background density and 
gas pressure are uniform and labelled as po and po respec- 
tively. However, we note that a spatial variation in po can cause 
phase mixing ([Heyvaerts & Priest 1983||De Moortel et al. 19991 
IHood et al. imW . 
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2.2. Coordinate system and non-dimensionalisation 

The linearised velocity, Vi, is split into components parallel and 
perpendicular to the equilibrium magnetic field. Thus, 



vi = V« 



Bn 



where Aq - -Bxz/a is the y-component of the vector magnetic 
potential. The terms in brackets are unit vectors. Splitting the ve- 
locity into components parallel and perpendicular to the equilib- 
rium magnetic field makes MHD mode detection and interpre- 
tation easier For example, in a low /5 plasma, the slow wave is 
guided by the magnetic field and has a velocity component that 
is mainly field-aligned. This makes perfect sense when yS «: 1 
but its usefulness is less clear when /? » 1 near the null point. 
To aid the numerical calculation, our primary variables are con- 
sidered to be Vx - VBo ■ BoVj. and vy = VBo ■ BoVy. 
Taking v, = means we do not consider the Alfven wave here as 
the description by 'McLaughlin & Hood (2004) is still valid (the 
plasma pressure plays no role in its propagation). 
We non-dimensionalise all variables by setting Vi - vovf, Vj^ - 
= V()Bv*, B() = BB*, Bi = BBt, x = ax*, z 

V = iV*, t 



Vofivl,V|| 



Pi = Pop 



0' 



aBAl and t] 



az , 

rio, where we 



let * denote a dimensionless quantity and vo, B, a,pQ, to and t/o 
are constants with the dimensions of the variable they are scal- 
ing. We then set —^= = Vo and Vn - — (i.e. we measure our 

speed in units of vo, which can be thought of as a constant back- 
ground Alfven speed). We also set ^ = 7?,^/, where R„ is the 

magnetic Reynolds number, and set ySo - where /3o is the 
plasma y8 at a distance unity from the origin (see Section l24b . 
Thus, we generate the non-dimensionalised versions of equa- 
tions (|2]i, ^ and (|5]) and under these scalings, t* - I (for ex- 
ample) refers to f = fo = i.e. the time taken to travel a dis- 
tance a at the reference background Alfven speed. For the rest 
of this paper, we drop the star indices; the fact that they are now 
non-dimensionalised is understood. 



2.3. Linearised equations 
The linearised equations are: 



dv± 
~di' 
5v|| 

dt 
db, 
~dt 
db, 
'df 
dpi 

dt 



1, Jdb, 
v,(.,z)^^ 



db, 
dpi 



2 \ dx dz 



_Pq Idp^ 

2 \ dx dz 



5v^ J_ 

dz Rm 

dv^ J_ 

dx R,„ 



d^b, d\. 



dx^ 

d^k 



dz^ 

d^b 



— i + ^Llll 

dx^ dz} 



-J 



x^ +z} 



5v 

5Vj_ 



dx dz 



-2 



x2 + z2 ''ll 



(z—+x 
\ dx dz 



4xz 
x^ + z^ 



(6) 



where the dimensionless Alfven speed, va (x, z) - Vx^Tz^. 
Note that by taking ySo = and ignoring resistiv- 
ity, (Rm oo), we recover the equations discussed in 
[McLaughlin & Hood (2004)1 



2.4. Piasma/3 

The parameter of key importance in equations (|6]l is ySg. The 
plasma /3 is defined as the ratio of the thermal plasma pressure to 
the magnetic pressure. In most parts of the corona, the plasma /3 
is much less than unity but, near null points, the magnetic field 
strength is small (and is zero at the null) and the plasma /? be- 
comes large. There is also coupling between the perpendicular 
and parallel velocity components when /So +- and this cou- 
pling is most effective where the sound speed and the Alfven 



velocity are comparable in magnitude. |Bogdan ef a/. (2003) call 
this zone the magnetic canopy or the /? 1 layer. 
The plasma fi varies throughout the whole region and 



2/Ufo 1 



y8 = 



A) 



(7) 



where - x^ + z^. Therefore, the fi - I layer occurs at a ra- 
dius r - "v^, where the gas pressure is equal to the magnetic 
pressure. However, in a uniform plasma, it is not the P - I 
layer that is most important, but instead it is where the sound 
speed is equal to the Alfven speed, i.e. Cs - va. Recalling that 

= (which impl 



IPo in non- 



lies Cj - 

dimensional units), we see that for the present non-uniform equi- 
librium Cj = at a radius r = yj^^ and it is when the incoming 

wave passes through this layer that the mode coupling occurs. Of 
course, the difference between the P - I layer at r = and 

the = VA layer at r = ^ is very small, and, hence, it is 
easier to refer to this as the y6 » 1 layer. 

Finally, we note that the basic fast and slow wave speeds for this 
equilibrium are given in dimensionless form as 



2 ^ r2y^Q/2 

3. MHD wave propagation with ySo - 



(8) 



0.25 



In this section, the linearised MHD equations, namely equa- 
tions (|6]l, are solved numerically using a two-step Lax-Wendroff 
scheme. In dimensionless units, the numerical domain is -4 < 
X < 4 and -4 < z < 4 and initially we consider a single wave 
pulse coming in from the top boundary due to a disturbance in 
the perpendicular component of the velocity. For the single wave 
pulse, the boundary conditions are chosen such that: 



(9) 



Vx(x,4) 




sin ojt for < f < 


n 


V||(x,4) 




otherwise 






= 


dv. 


dv± 


z=-4 = , 


-:^l.v=-4 
dx 


dx 


IT 






dvu 






-3— Lt=-4 

dx 


= 


, ^L=4 = , 
dx 


dz 


=-4 = 0. 



Tests show that the central behaviour is largely unaffected by 
the choices of side and bottom boundary conditions. The other 
boundary conditions on the perturbed magnetic field follow from 
the remaining equations and the solenoidal condition. 
The location of the upper boundary is not important and it is 
the value of /3q that determines the distance the fast wave trav- 
els before encountering the mode conversion region. Initially, 
we set A) - 0.25 and R^ - 10^. We also take oj - 2n, as in 
McLaughlin & Hood (2004)| The results for v^ can be seen in 
Figures |2] and[3Jand for vy in Figure |4] 
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3. 1. Initial disturbance in the perpendicular component 

In the zero j6 limit (see McLaughlin & Hood 2004), the wave 
pulse given above is a fast MHD disturbance, but for the present 
problem this is no longer true. In the system presented here, the 
fast wave is no longer purely described by a perpendicular com- 
ponent: it is predominately given by the perpendicular compo- 
nent but also has a smaller parallel component. Thus, the single 
wave pulse described above is a disturbance in tlie perpendicu- 
lar component, as opposed to simply a fast MHD disturbance. 
This difference in definition will be important when we come to 
interpret the different wave types in SectionH] In addition, it can 
be seen that the governing equations (|6]l that v^ acts as a driver 
for V||. Since the value of the plasma p is small on the upper 
boundary, v^ initially behaves in the same manner as the zero 
jS case. However, the parallel component of the velocity is now 
driven by the perpendicular component. 

Thus,the general solution for V|| will consist of two parts; a com- 
plementary function, corresponding to a slow mode disturbance 
and a particular integral, due to the fast mode driver term. These 
two parts to the vy wave are clearly seen in the simulations. 
However, it is clearer if we first discuss the behaviour of the 
Vj. component before moving on to the parallel component. 



3.2. Behaviour ol ^ ±_ 

We find that the v^ disturbance travels towards the neighbour- 
hood of the null point and begins to wrap around it. This is 
due to refraction caused by the spatially varying Alfven speed, 
va(x,z), as noted by |Nakariakov & Roberts (1995)| and seen in 
McLaughlin & Hood (2004) (further references are cited in the 
latter). Once the wave reaches the = va layer (denoted by the 
black circle around the null in Figure |2]l the nature of the wave 
changes; part of the wave now appears to spread out along the 
field lines (this part also moves slower than the rest of the wave). 
Meanwhile, the majority of the wave outside the circle contin- 
ues to refract about the null. This wrapping eff'ect repeats, and 
each time part of the wave crosses the - va layer, part of it 
is transformed to the field-guided wave. Figure [3] shows a blow- 
up of subfigures {d) to (/) during the time of the first crossing. 
In this figure, we can clearly see the splitting of the wave as it 
nears and then crosses the c, - va layer. Thus, we can see that 
there are three main types of wave behaviour in the perpendic- 
ular component; outside the = va layer we predominatrely 
see the refraction effect, whereas inside we see both fast and 
slow wave behaviour; the latter which appears as a field-guided 
wave. We shall give a fuller description of this phenomenon in 
Section|6]and give an interpretation of these different wave types 
in Section]?] 



lytic solution (see Appendix A). There is also a complementary 
function part that satisfies the homogeneous equation 



3.3. Behaviour of \\\ 

The behaviour of vy can be seen in Figure ]4] We see that it is an 
antisymmetric wave about the z axis, and that it is focusing in 
towards the null. This pulse travels at the same speed as the per- 
pendicular component wave, namely the fast wave speed. This is 
not surprising since, as mentioned previously, Vj^ drives vy. Thus, 
mathematically this leading wave can be thought of as the par- 
ticular solution to an (inhomogeneous) equation. Note that the 
range of the colour scale indicates that the maximum amplitude 
of V|| is 0.2 and this is approximately the same size as jSq. We use 
this relatively weak coupling to construct an approximate ana- 



C?2 



(Bo ■ V)2 



CF 



dt^ 2 

This has a solution that involves a function of the variable 



(10) 




The complementary function part can be seen in Figure ]5] This 
second part of the solution lags behind the particular solution and 
propagates at the slow wave speed. In a low p plasma this is ap- 
proximately the sound speed and hence is proportional to yf^Q. 
From the lower subfigures of Figure ]5] we see that the comple- 
mentary function does indeed travel with a speed proportional 
to yfPa, and has an amplitude that varies as fi^^. The comple- 
mentary function has a discontinuous shape (as is also shown in 
Appendix A) and its amplitude is small, 0(J3q^), compared to the 
particular integral, 0(J3o), and is thus difficult to see in Figure ]4] 
There are also several other noteworthy aspects to Figure ]4] 
Firstly, it is obviously more complicated than the Vj^ wave. 
Secondly, this component has eight lobe-like structures by the 
time f = 2.0, when the leading edge of the wave is about to cross 
the region where c, - va. The lobes alternate between positive 
and negative values. This complicated lobe-like structures can be 
explained by considering the same system but driven by a circu- 
lar boundary condition. This is explained in Appendix A, where 
the main conclusions is that the equilibrium magnetic config- 
uration naturally leads to a sin 40 dependence in vy. This ex- 
plains the complicated, lobe-like structuring, i.e. the lobes come 
from the equilibrium magnetic configuration (but the wave is no 
longer circular so the sin 46* is less apparent). 

4. Interpretation 

Roberts (1985) showed that VAifven, Vsiow and Vfast form an or- 
thogonal basis of vectors for the linearised MHD equations in 
a uniform plasma. In this paper, we do not consider the Alfven 
wave and so our 2D system can be described in terms of the vec- 
tors Vfast and Vsiow Due to the form of the equilibrium magnetic 
field near the null point, we choose to work in the directions par- 
allel and perpendicular to the magnetic field and thus we may 
represent these two vectors in terms of Vfast and Vsiow, i-S- 

y± = ^Vfast + SVslow , V|| = CVfast + ©Vsiow 

where S, C and D are unknown functions that depend 
upon the magnetic geometry (and possibly the plasma /3). 
Equivalently, we may represent Vfast and Vsiow in terms of and 
V||. This representation is only possible because both Vfast & Vsiow 
and & vy are linearly independent vectors. Note that in a low 
yS region or when the plasma is nearly uniform, S will be very 
small. Thus in these regions, we interpret as predominately a 
fast wave (i.e. Vj. ~ ^Vfast). 

However, we must be careful. The concepts of fast and slow 
waves were derived for a unidirectional magnetic field (see 
Edwin & Roberts 1983) and so these ideas may not carry over 
to more complicated geometries quite as simply as claimed here. 
Nonetheless, we shall continue to use terms such as fast and slow 
wave in our interpretation of the waves in this paper (i.e. we shall 
use the terminology and intuition gained from the unidirectional 
magnetic field model). Thus, we choose to interpret tlie waves 
seen in our perpendicular and parallel velocities using the ter- 
minology of fast and slow waves. 
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Fig. 2. Contours of for numerical simulation for a wave sent in from the upper boundary for -4 < x < 4 and /3o - 0.25 and its 
resultant propagation at times (a) f=0.33, (b) f=1.0, (c) f=2.0, (d) f=2.33, (e) f=2.67, (/) f=3.0, (g) t =3.33, (/z) f=3.67 and (/) f=4.0, 
labelling from top left to bottom right. The black circle indicates the position of the Cs = va layer The cross denotes the null point 
in the magnetic configuration. 
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Fig. 3. Blow-up subfigures of from Figure|2]at times (a) f=2.33, (b) t-2.61 and (c) f=3.0, labelling left to right. 



5. Mode conversion across the Cs = va layer 

The Cs = Va layer is of critical importance to our system. When 
the Alfven speed and sound speed are dissimilar, there is negligi- 
ble coupling between the fast and slow magnetoacoustic waves. 
However, near the c., - va layer, the two waves can resonantly 
interact with each other and strong mode coupling can occur. 
Since /3 — -j^^ (where /3q is a constant of our choosing), our 
system consists of a region of low /3 plasma outside the c, - va 
layer and a region of high /3 plasma within (see Figure |6]l. Fast 
and slow waves have differing properties depending on if they 
are in a high or low /3 environment (see e.g. Bogdan et al. 2003). 
See Table 1 for a summary of these properties. 
Thus, in a low /3 region, the velocity vector of the fast wave is 
predominately perpendicular to the equilibrium magnetic field. 




-3 -2 -1 



Fig. 6. Regions of high and low jS in our equilibrium magnetic 
field, where /? - -^rp'- The black circle indicates the position of 
the Cs — Va layer and the cross denotes the null point. 
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Fig. 4. Contours of vy for numerical simulation for a fast wave sent in from upper boundary for -4 < x < 4 and /3o = 0.25 and its 
resultant propagation at times (a) t=033, (b) t=1.0, (c) f=2.0, (d) f=2.33, (e) t=2.61, (/) f=3.0, (g) t =3.33, (h) f=3.67 and (/) f=4.0, 
labelling from top left to bottom right. The black circle indicates the position of the = va layer The cross denotes the null point 
in the magnetic configuration. 



although the wave can propagate almost isotropically. The slow 
wave has a velocity component along the field and propagates 
along B(). In the high /? region, the fast wave again propagates 
isotropically but this time the velocity vector is parallel to the 
wavenumber vector Hence, the fast wave now has components 
of velocity that are perpendicular and parallel to the field. The 
slow wave is again guided by the field but the velocity is perpen- 
dicular to BoQ 

In our numerical simulations, a straight wave pulse in the per- 
pendicular component is sent in from the low /3 region at the 
top boundary. This disturbance is a low /3 fast wave. At some 
point the wave will cross the c, = layer and enter the high 
j3 environment. Because of the above properties, the incoming 
fast wave is initially perpendicular to Bo but, on crossing the 
Cj = Va layer, part of this perpendicular component remains as 
a fast wave and part is converted into the slow wave. The first 
point of contact occurs along x = and at later times for in- 
creasing values of x^. Thus, we have a low p wave approaching 
the layer, coupling and mixing inside the layer and emerging as 
a mixture of high /3 fast and slow waves. Later on, the fast wave 
leaves the high /3 region and re-enters the low /3 region. Mode 
conversion can occur at this crossing as well. 



' Note that the high /? slow wave shows up in the perpendicular com- 
ponent because it is a transverse wave (v • k = 0) and not because it is 
propagating across the field; the slow wave cannot propagate across the 
field! 



6. MHD wave propagation with j3o = 2.25 

In this section, the linearised equations (|6]l are again solved nu- 
merically but for fSo - 2.25 so that the c, = va layer is now at a 
larger radius than above. All other boundary conditions remain 
the same. The behaviour of the perpendicular component of the 
magnetoacoustic wave when it is inside the c, = va layer can be 
seen in Figure|7] Again, we find that the linear, magnetoacoustic 
wave travels towards the neighbourhood of the null point and be- 
gins to wrap around it (as it did forySo - 0.25). In this case, since 
the yS ~ 1 is now reached earlier, the wave has had less time to 
refract. Using our terminology and interpretation from Sections 
|4]and|5] we identify this as the low j3 fast wave. However, when 
the wave passes through the C j = va layer, a secondary wave be- 
comes apparent; we identify this as the high /5 slow wave. Thus, 
the low /3 fast wave has transformed into two high /3 waves (the 
split can be followed closely in Figure|7]i. 
The propagation now proceeds in three ways: 

- The generated high /3 slow wave spreads out along the field 
lines and accumulates along the separatrices. This slow wave 
cannot cross the separatrices. 

- The high p fast wave continues to refract and some of the fast 
wave (located very close to the null point) passes through the 
null. This efifect was not seen in the /5 - model. The fast 
wave can now pass through the null because there is now a 
non-zero fast mode speed, c fast, at the origin. Thus, the high 
P fast wave has passed (slowly) through the null. This cross- 
ing still creates a large accumulation of current near the null. 
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Complementai'y function 




Fig. 5. Shaded surface of vy for numerical simulation for a fast wave (pulse) sent in from upper boundary for -4 < x < 4 and 
/3o = 0.25 at time f = 1.5 (top left) and f = 3.5 (top right). The bottom left figure shows the position reached by the complementary 
function along x = 1.5 after f = 1.5 (the gradient of the straight line based on the first and last points is -0.98. The bottom right 
figure shows the amplitude of the complementary function as a function of ySo. The amplitude is proportional to although /3^-^ 
is to be expected from the boundary conditions. 



However, unlike the /3 - case, it remains finite. The high /3 
fast wave continues to propagate downwards and leaves the 
Cs = va layer, converting to a low /3 fast wave and spreading 
out (low j6 fast wave propagates roughly isotropically). 
- Finally, the rest of the fast wave located away from the null 
(low j6 fast wave) is not greatly affected by the non-zero 
sound speed (as » '^sound-' continues to refract 

around the null. In fact, as the wings of the low /3 fast wave 
wrap around below the null point, they encounter the part 
of the fast wave that has travelled through the null and is 
leaving the c, - va layer . This results in a complicated in- 
terference pattern. Nonetheless, it appears that the two fast 
waves passes through each other without influencing each 
other (due to the linear nature of the system). A full non- 
linear treatment of the equations may reveal a different be- 
haviour. 

Thus, we have a different behaviour than that seen for the j6o - 
0.25 case: With a large /3q, the perpendicular component of the 
velocity has passed through the null (interpreted as a high /3 fast 
wave crossing the null), whereas with a small ySo, we did not 
see any significant wave pass through the null by the end of the 
simulation. 

So what is the relationship between the fio = 2.25 and /3o = 0.25 
cases; are they not just scaled versions of each other? This will 
be discussed in the conclusions. 

6.1. Quantifying mode conversion 

At present, there does not exist a robust set of rules connecting 
low and high /3 waves across the Cj - va layer or across the y6 = 1 
layer dCarlsson & Stein 19991 ). Here we specifically investigate 
what happens when a low /3 fast wave crosses the Cj = va layer 
and becomes part high fast wave and part high /3 slow wave 
(or alternatively what happens when it crosses the a: 1 layer a 
little time earlier). An obvious question to ask is how much of the 



incident wave is converted to high /3 slow wave. This is difficult 
to quantify in our /3o = 0.25 simulation, as the high /3 fast and 
slow waves never really separate enough for us to measure them 
individually. However, the waves do completely separate for the 
/3() - 2.25 simulation, for example in subfigure 7(g) the slow 
wave is clearly separated from the other wave types (i.e. the high 
and low /? fast waves). Thus, we can measure the amount of this 
slow wave part by integrating v^ over the area of the slow wave 
(to get the volume or, since the equilibrium density is uniform, 
the total momentum of the slow wave) and compare it to the 
original incident wave pulse (given by the boundary conditions; 
equations (|9]l). This shows that 9.2% of the initial low fast 
wave disturbance is converted to high /3 slow wave (for the given 
boundary conditions with /?o - 2.25). We can also work out how 
the proportion of initial disturbance converted to high /3 slow 
wave changes with /3o. This is shown in Figure |8] where we can 
see that the fraction of original disturbance converted to high /? 
slow wave is proportional to f3o. 




o - 
5 r 



^2 4 6 8 10 

Fig. 8. Fraction of original disturbance converted to high fi slow 
wave against parameter y6o, (gradient of straight line is approxi- 
mately 2.0). Note that for ySo < 2.0, the two wave types do not 
separate enough for us to measure them individually. 
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Fig. 7. Contours of Vj. for numerical simulation for a fast wave sent in from upper boundary for -4 < x < 4 and Pq - 2.25 and its 
resultant propagation at times {a) t^lJ, (b) ?=1.8, (c) f=1.9, (d) t = 2.0, (e) t^l.l, (/) t = 2.2, (g) f=2.3, {h) t=2A and (/) f=2.5, 
labelling from top left to bottom right. The white circle indicates the position of the = Vyi layer and the cross denotes the null 
point in the magnetic configuration. 



7. Conclusions 

This paper extends the investigations of McLaughlin & Hood 
(2004) into the nature of magnetoacoustic waves in the neigh- 
bourhood of a null point. We have conducted two experiments, 
varying the choice of /3o in each, where the choice dictated the 
location of the Cj = va layer. 

In the first experiment (Section |3]l, we set j3o - 0.25. We find 
that the linear magnetoacoustic wave travels towards the neigh- 
bourhood of the null point and begins to wrap around it. This oc- 
curs due to the spatially varying Alfven speed, ix,z)- +z^, 
and has been observed before (Nakariakov & Roberts 1995; 
McLaughlin & Hood 2004). In this experiment, we drive a wave 
pulse in the perpendicular velocity component and we identify 
this pulse as a low /3 fast wave. However, once the wave reaches 
and crosses the c , - va layer, part of the wave is transformed into 
a high p slow wave and the majority into the high p fast wave. 
The slow wave part spreads out along the fieldlines. Meanwhile, 
the rest of the (low P) fast wave continues to wrap around the 
origin. This refraction effect continues again and again, but each 
time part of the fast wave is converted to a slow wave as it crosses 
the Cs = Va layer The majority of the current build-up occurs 
very close to the null. 

In our second set of simulations (Section |6]l, we set ySo - 2.25, 
thusplacing the Cj - va layer at a larger radius than before (since 
in this magnetic geometry, the c, - va layer occurs at a radius 

r - a/^). In this experiment, we find that the magnetoacous- 



tic wave travels towards the neighbourhood of the null point and 
begins to wrap around it (as before). However, when the wave 
crosses the Cj = va layer, a secondary wave (high p slow wave) 
becomes apparent (we find that the fraction of incident wave 
converted to slow wave is proportional to P). The propagation 
now proceeds in three ways. Firstly, the generated slow wave 
spreads out along the fieldlines and accumulates along the sep- 
aratrices. Secondly, the remaining part of the fast wave inside 
the Cs = Va layer continues to refract and some of it (located 
close to the null point) passes through the origin. We identify 
this part as a high p fast wave. The high p fast wave can pass 
through the origin because, although va{0, 0) - 0, there is now a 
non-zero sound speed there (and c^^^^ ~ ^1 + '^sound^- ^^^^ pass- 
ing causes a large current accumulation near the origin. After it 
has crossed the null, the high p fast wave continues downwards 
and leaves the Cj - va layer. As it emerges, it becomes a low 
P fast wave and spreads out (since the low p fast wave propa- 
gates almost isotropically). Finally, the fast wave located away 
from the null and c, = va layer (the 'wings' of the low p wave) 
are not affected by the non-zero sound speed (as » '^sound^ 
and so for them the refraction eff'ect dominates. In fact, as these 
wings wrap around below the null point, they encounter the high 
P fast wave as it is emerging from the Cj = va layer This re- 
sults in a complicated interference pattern, but it appears that 
the two waves passes through each other without any lasting ef- 
fect on each other (due to the linear nature of the system). The 
Pq - 2.25 numerical simulations also showed good agreement 
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with a WKB approximation, until the breakdown point when the 
neighbouring rays cross each other 

The part of the fast wave that goes through the null effectively 
escapes the refraction effect of the null point. It seems as if the 
smaller the value of ySo used, the slower the high fi fast wave 
can cross the null. This can be understood since the high fi fast 
wave travels close to the sound speed, which is related to jSq 

(recall from section l24l that c* = sj^Pa)- Thus, the behaviour 

of the j3\^ - 2.25 and fSo - 0.25 cases can be understood; in 
the j8o = 0.25 case, the high p fast wave can pass through the 
null, but it travels at such a slow speed (as c , varies as -v^) that 
the refraction effect in the low fi region dominates. Hence, as the 
Cj = va_ layer gets closer to the null, less and less of the fast wave 
can pass through the null in a given time, and the parts that can 
travel through do so at a slower and slower speed. Also, these 
parts of the wave that do pass through tend to be be swamped by 
the wings of the rest of the wave wrapping round below the null 
(repeatedly). 

This explains the relationship between our y6o = 0.25 and 2.25 
investigations. These experiments show scaled versions of each 
other, but with the larger value of ySo it is clearer to see that 
there are actually two competing phenomena; a refraction ef- 
fect caused by the varying Alfven speed and a non-zero sound 
speed at the null which allows the fast wave to pass through. It 
is the value of Pq that dictates which effect dominates. 
Thus, two extremes occur. The first occurs when Po — > 0; in 
which case the refraction effect infinitely dominates over the 
other effects (and we recover the results of McLaughlin & Hood 
2004) and the second when /Jq — » oo and the system becomes 
hydrodynamic. In this case, the fast wave reduces to an acoustic 
wave and so completely passes through the null (in effect it does 
not even see the magnetic field, since <K c,). Thus, we can 
understand the whole spectrum of values of the parameter ySo. 
From this work, it has been seen that a warm plasma introduces 
many new effects not seen in the cold plasma limit, most notably 
the introduction of fast and slow waves to the system. It also 
appears that the choice of Pq is of critical importance, since the 
two experiments yield different results. However, the key choice 
here is not in picking Pq, since this just determines where the 
Cj = va layer will occur. The choice of Pq is equivalent to choos- 
ing where to set the boundaries of our box, and so the critical 
parameter in our system is in choosing the distance between the 
initial disturbance and the c^ — va layer. This is because these 
experiments do not simply show scaled versions of each other, 
they show two competing phenomena; a refraction effect caused 
by the varying Alfven speed, and a non-zero sound speed at the 
null which allows the fast wave to pass through. 
If the c , - Va layer is close to the null (small choice of Pq and 
hence there is not much coupling to the parallel velocity and 
pressure terms) then the refraction effect will dominate, result- 
ing in a scenario similar to Po - 0.25, the waves remain trapped 
near the null point and heating will occur close to the null. The 
extreme case is Pq — > 0, in which case the refraction effect in- 
finitely dominates over the other effects and we recover the re- 
sults of IMcLaughhn & Hood (2004) 



However, if the Cj = va layer is far from the null, a portion of 
the fast wave will be able to pass through the origin and escape 
the system (i.e. will not deposit its energy near the null point). 
In this case, the remaining wave energy of the generated slow 
waves is dissipated along the separatrices. So heating will occur 
in both systems, although the nature will be different. 
This is all very interesting for its mathematical sake, but how 
does the competition of these two effects manifest itself in the 



corona? The plasma p parameter is defined as the ratio of the 
thermal plasma pressure to the magnetic pressure. In most parts 
of the corona, /? <k 1 and hence the pressure gradients in the 
plasma can be neglected. Values of p - 0.01 are often quoted 
(e.g. Priest 1982). However, near null points the plasma p can be- 
come very large, so the (true) plasma p varies through the whole 
region. However as seen above, it is the distance between the 
initial (planar) pulse and the c, - va layer that is of critical 
importance. We believe that coronal disturbances will propagate 
for some distance before they encounter a (coronal) null point. 
Hence, the fast wave part of the disturbance will feel the refrac- 
tion effect of the null and begin to refract around it. By the time 
the c, - Va layer is reached or the sound speed becomes impor- 
tant, the fast wave disturbance will be (almost) circular in nature 
and will have thinned and dissipation will be extracting the en- 
ergy from the wave. It is true that some of the wave may pass 
through the null or be converted into high p slow waves, but it is 
likely that the majority of the wave energy will accumulate close 
to the null, causing large current accumulation and heating there. 

Appendix A: Lobe generation 

We can gain insight into the p system using a small Pq expan- 
sion of the linearised ideal MHD equations. In polar coordinates, 
our magnetic field is 

Bo = -r cos 261 f H- r sin 20 

We now assume v^ and vy can be expanded in powers of Pq such 
that; 

= Vi 1 (r, f) + Pi ^±2(r,0,t), V|| = Pq Vy 1 (r, 0, t) . 
Substituting these forms into equations (|6|l gives: 



5f2 

dt 
d_ 
It 



-OifiQ) 



V|li = 



d_l d_ 

dr\ dr ^' 

-]^{^Q ■V)p,+ OifiQ) , 



-yV ■ 



sin 20 



cos 20 



+ 0(fiQ) 



- y sin 



(11) 



Substituting into the equation for ^V||j and ignoring terms 
of order gives: 

CJ2 



dt^ 



J ( d 
— — < -rcos20sin20-— r 
2 I dr 



dr\ r^ I 



+ r sin 20 cos 20 



y 



d 1 



- sin 40 r 

4 drr 



d /v^n 



^ sin 40 T(r,t) 



(12) 



where "Fir, t) is a known function that depends only on r and 
t. This function is determined solely by the behaviour of v^ i . 
Although the shape of the v^i pulse may be continuous, the 
derivatives in (fT2l l mean that vy may have discontinuous lead- 
ing and trailing edges. Thus, we can see that, with our choice of 
magnetic null point and by driving v^, the form of vy naturally 
develops a sin 40 dependence. Assume an initial circular pulse 



sin;r(r - 2.5) 2.5 < r < 3.5 







elsewhere 
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Fig. 9. (Left) Contour of | sin 40^ [(s - n^r^) sin n(r - 2.5) -Snr cos 7r(r - 2.5)] for 2.5 < r < 3.5. (Right) Contours of the numeri- 
cal simulation of pi for a fast wave pulse initially located about a radius Vx^ + - 3 and its resultant propagation after time t-l.O. 
The black circle indicates the position of the c , = va layer and the cross denotes the nuU point in the magnetic configuration, pi has 
an alternating form, where orange represents pi > and blue pi < 0. 



The right hand side of ( fT2l i is 
|sin46i^ [(8 -7rV2)sin7r(r-2.5) -57rrcos7r(r-2.5)] and a 
contour plot of this is shown in Figure |9] However, the contour 
shown is not a true integral of such an initial condition (i.e. 
is not a solution of (fT2] i) but it does clearly demonstrate the 
sin 40 dependence. From equation (fTTT i. we can also see that 
p\ develops a sin 20 behaviour. This behaviour can be seen on 
the right-hand side of Figure |9] which shows the contours of p\ 
from a numerical simulation. 



Appendix B: Analytical work 

In order to gain more insight into the numerical simulations, we 
use a geometrical optics WKB solution to obtain an approximate 
analytical solution. By combining equations Q and dU, we 
can form: 



di- 



= 2Z£v(V-v) + {Vx[Vx(vxBo)])x^ 

Po 



(13) 



We substitute v - ae"^^-^'^^ ■ e^"^' into ([13]) where a is a constant 
amplitude vector in the geometrical optics approximation. The 
physical optics approximations would allow the amplitude vec- 
tors to vary with position as well as the wavenumber (p. Thus, we 
follow the ray paths. We now take the dot product with Bo and 
V0 to get the two velocity components: 



0) 



-^(B„-V0) 



vBo 
V ■ V0 



These two coupled equations must have zero determinant to pre- 
vent a trivial solution, and so by evaluating the determinant, sub- 



stituting Cj = ^ and Vq = and making the WKB approxi- 
mation such that (jj ~ (p^ \ leads to a first order equation of the 
form: 

+ c]vl{p^+q'-)(xp-zqf 
^ 'F(x,z,<p,p,q) = 

= l[co^-co'[cl+yl{x'+Z%p' + q') 

+ C?v2(/ + ^2)(x;,-z^)2] 



where p - j^. and q - '-^ and ;F is a non-linear PDE. Also note 
that the Alfven speed - Vq {x^ + z^). This PDE can also be 
written as 

2cJ^[c]+.l{x^+z')]{p'+q') 
± ^^ (,p^ + q^f \cl + \l (x^ + z^)f - 4 c^v^ (p^ + q^) (xp - zqf 

This equation is reminiscent of the dispersion relation for mag- 
netoacoustic waves (e.g. Roberts 1985). This equation contains 
information about the two wave types. The method to solve this 
equation is to assume it to be of the form (u^-(i?y^^)((jip--(jL>^^^^ - 
0. Considering the fast wave (so ap- + w^j^^), we can apply the 
method of characteristics to generate the equations: 



ds 
dp 

Ys 
dq 

ds 
dx 

ds 
dz 
ds 



2co^ 

-Ax - [CA^x - IpABc^vl) ID 

-Az - [CA^z + 2qABc]v^^ /D 
Cp + [AC^p - 2pB^c]wl - 2cIwIxAB) jD 
Cq + {AC^q - 2qB^c]vl + 2c]vIzAB) ID 



where A 



p2 + ^2, 



B - xp - zq, C 



cl + vl 



x^ + z^), D 



4 



— AAB^c]\-^ and to is the frequency of our wave and s is 

some parameter along the characteristic. These five ODEs were 
solved numerically using a fourth-order Runge-Kutta method. 
Contours of constant (p can be thought of as defining the posi- 
tions of the edges of the wave pulse, i.e. with correct choices of 
s, the WKB solution represents the front, middle and back edges 
of the wave. 

We see that the agreement between the numerical simulation and 
the analytical approximation is quite good (see Figure [TOl i. In 
Figure[TT] we see the main features of the numerical simulation; 
the wave refracts as it approaches the null and that part of the 
wave passes through the null (at a slower speed than the 'wings' 
outside the Cj - va layer). This central part (that passed through 
the null) then emerges from the = va layer and spreads out 
isotropically. Meanwhile, the wings continue to refract around 
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Fig. 10. Comparison of from numerical simulation and analytical solution for a fast wave sent in from upper boundary for 
-4 < ;c < 4 and ySo = 2.25 and its resultant propagation at times (a) f=0.33, (b) r=1.0, (c) f=2.0, (d) f=2.33, (e) f=2.67, (/) f=3.0, 
(g) t =3.33, (h) t=3.61 and (/) f=4.0, labelling from top left to bottom right. The lines represent the front, middle and back edges of 
the wave. The circle indicates the position of the c., = va layer and the cross denotes the null point in the magnetic configuration. 



the null. Also, as the wave crosses the c, — va layer, the wave- 
front overlaps with itself (forming small triangular shapes, see 
FigurefT2l): this may explain why we see the perpendicular wave 
form a sharp edge as it crosses the - va layer. However, 
the crossing of neighbouring rays frequently indicates the break- 
down of the WKB solution, as happens at caustics. This is almost 
certainly where the conversion of part of the fast wave into a 
slow wave occurs and this process requires a more detailed study 
of the WKB equations. Instead of the triangular shape propagat- 
ing at the fast speed it should propagate at the slow speed. 
The rays of the WKB solution can be seen in Figure[T2l The left 
graph shows the smalltriangular shape formed as the wavefront 
overlaps with itself when it crosses the = va layer, indicating 
that this is not just a discontinuity in the wavefront (i.e. since it 
can be resolved). The central graph shows the rays for starting 
points of X = 1,2, 2.5, 3 and 4 along z - 4. We see that for a 
starting point of x = 3 and x = 4 (z = 4), the ray is deflected sub- 
stantially by the null point, whereas the deflection is less severe 
for 2 and 1. A starting point of x = 2.5 seems to be the critical 
starting point that determines if ^ changes sign (at least in the 
range -4 < x < 4). We also see that x = is not deflected at all. 
The right hand graph shows the rays for for x e [0, 4] in inter- 
vals of 0.1. The dashed line represents x = 2.5; for x > 2.5, the 
ray path is deflected so much that ^ changes sign. For x < 2.5, 
the rays are deflected but escape from the opposite corner of the 
system. This nature reflects the numerical solution; part of the 
wave passes through the null and part refracts around it. 



However, the analytical WKB approximation does not give a full 
description of the numerical simulation. Firstly, the WKB ap- 
proximation breaks down when neighbouring rays cross and this 
requires a more detailed study to demonstrate exactly how the 
fast and slow waves interact. Secondly, the WKB approximation 
can also be applied to the a/ - (^^^^^ = equation but this does 
not greatly add to our understanding of the system. 
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Fast Wave 


Slow Wave 


High^e 


Behaves like isotropic sound wave 
(speed Cj) 


vJ LIlUvl^Ll tXlKJll^ '-'0 

Transverse wave propagating at speed Va 
Dominant velocity component is v^ 


Low>8 


Propagates roughly isotropically 
(speed Va) 
Dominant velocity component is Vj^ 


Guided along Bo 
Longitudinal wave propagating at speed Ct 
Dominant velocity component is Vy 



Table 1. Properties of fast and slow waves depending on their environment. 



